Observation of Phase Defects in Quasi-2D Bose-Einstein Condensates 
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We have observed phase defects in quasi-2D Bose-Einstein condensates close to the condensation 
temperature. Either a single or several equally spaced condensates are produced by selectively 
evaporating the sites of a ID optical lattice. When several clouds are released from the lattice and 
allowed to overlap, dislocation lines in the interference patterns reveal nontrivial phase defects. 
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Low dimensional bosonic systems have very different 
coherence properties than their three dimensional (3D) 
counterparts. In a spatially uniform one dimensional 
(ID) system, a Bose-Einstein condensate (BEC) cannot 
exist even at zero temperature. In two dimensions (2D) 
a BEC exists at zero temperature, but phase fluctua- 
tions destroy the long range order at any finite temper- 
ature. At low temperatures the system is superfluid and 
the phase fluctuations can be described as bound vortex- 
antivortex pairs. At the Kosterlitz-Thouless (KT) tran- 
sition temperature the unbinding of the pairs be- 
comes favorable and the system enters the normal state. 

In recent years, great efforts have been made to study 
the effects of reduced dimensionality in trapped atomic 
gases 0. In both ID and 2D, the density of states in a 
harmonic trap allows for Bose-Einstein condensation at 
finite temperature. In contrast to ID and elongated 3D 
systems [illlllili[lS[ni[Il[li,the coherence prop- 
erties of 2D atomic BECs have so far been explored only 
theoretically [T ^ [l^ ll^. In previous experiments, 
quasi-2D BECs [d llSl Il9l l20| or ultracold clouds 
were produced in specially designed "pancake" trapping 
potentials. The sites of a ID optical lattice usually also 
fulfill the criteria for 2D trapping jH El [3 H^; the 
difficulty in these systems is to suppress tunneling be- 
tween the sites, and to address or study them indepen- 
dently mil. 

In this Letter, we report the production of an array 
of individually addressable quasi-2D BECs. By selec- 
tively evaporating the atoms from the sites of a ID op- 
tical lattice, we can produce either a single or several 
equally spaced condensates. The distinct advantage of 
this approach is that it opens the possibility to study the 
phase structures in quasi-2D BECs interferometrically. 
We have observed interference patterns which clearly re- 
veal the presence of phase defects in condensates close to 
the ideal gas Bose-Einstein condensation temperature. 
We discuss the possible underlying phase configurations. 

Our experiments start with an almost pure ^''Rb con- 
densate with 4 X 10^ atoms in the F — ttif ^ 2 hy- 
perfine state, produced by radio-frequency (rf) evapo- 
ration in a cylindrically symmetric loffe-Pritchard (IP) 
magnetic trap. The trapping frequencies are ll!z/'2t^ = 
12 Hz axially, and lu±/2tt — 106 Hz radially, leading to 
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Figure 1: An array of individually addressable quasi-2D 
BECs. (a) A ID optical lattice splits a cigar-shaped 3D 
condensate into 15 - 30 independent quasi-2D BECs. (b) A 
magnetic field gradient along the lattice axis allows us to se- 
lectively address the sites by an rf field. We evaporate the 
atoms from all the sites except those within a frequency gap 
A. (c) Steps in the BEC atom number Aat as a function of 
A, corresponding to 0, 1 and 2 sites spared from evaporation. 
Each data point represents a single measurement. 

cigar-shaped condensates with a Thomas-Fermi length 
of 90 /im and a diameter of 10 ^m. 

After creation of the BEC we ramp up the periodic 
potential of a ID optical lattice, which splits the 3D 
condensate into an array of independent quasi-2D BECs 
(see Fig. ^a) and [l^l)- The lattice is superimposed on 
the magnetic trapping potential along the long axis (z) 
of the cigar. Two horizontal laser beams of wavelength 
A — 532 nm intersect at a small angle 9 to create a stand- 
ing wave with a period of d = X/[2sin{9/2)]. The blue- 
detuned laser light creates a repulsive potential for the 
atoms, which accumulate at the nodes of the standing 
wave, with the radial confinement being provided by the 
magnetic potential. Along z, the lattice potential has the 
shape V{z) = Vbcos2(7rz/d), with Vo/h w 50 kHz. 

For the work presented here we have used two lattice 
periods, di = 2.7 /xm and ^2 = 5.1 /xm. The respec- 
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tive oscillation frequencies along z are u}\/2it = 4.0 kHz 
and ujil^-n = 2.1kHz. At the end of the experimen- 
tal cycle (described below), the BEC atom numbers in 
the most populated, central sites are A^i « 10^ and 
7V2 « 2 X 10'*. We numerically solve the Gross-Pitaevskii 
equation to get the corresponding chemical potentials 
fii/h = 2.2 kHz and /X2//1 = 2.5 kHz, where the haJi^2/'2 
zero-point offset is suppressed in our definition of /i. In 
the smooth crossover from 3D to 2D, the condensates in 
the shorter period lattice are thus well in the 2D regime 
with fii/{huji) = 0.6, while for the clouds in the longer 
period lattice this ratio is 1.2. 

Since the radial trapping is purely magnetic, we can 
remove atoms from the lattice by rf induced spin-flips to 
untrapped Zeeman states. In order to address the lat- 
tice sites selectively, we apply a magnetic field gradient 
h' along z [23, HI] . This creates an energy gradient along 
the lattice direction, and splits the resonant frequencies 
for evaporation of atoms from two neighboring sites by 
= A*B&'t^i,2/(2/i), where /xb is the Bohr magneton 
(Fig.^b)). We use gradients up to 26 G/cm, correspond- 
ing to /S.vi = 5 kHz and Ai'2 = 9 kHz. These splittings 
are larger than the chemical potentials /Lti,2, and the rf 
Rabi frequency (« 2 kHz) . The lattice sites can thus be 
addressed individually. 

The experimental routine to produce an adjustable 
number of condensates starts with a slow, 200 ms ramp- 
up of the gradient b'. As illustrated in Fig.^b), we then 
evaporate the atoms from both ends of the cigar, sparing 
only the central sites within a variable rf frequency gap 
A. We perform this evaporation in 100 ms, switch off the 
rf field, and ramp b' back to zero in another 200 ms [i^ . 
During this time, some heating of the remaining clouds 
occurs, and they reach a temperature slightly below the 
condensation temperature, as we discuss in detail below. 

To verify that we can address the lattice sites individu- 
ally, we measure the total number of condensed atoms left 
in the trap as a function of A. An example of such a plot 
is shown in Fig. ^c) for d2 = 5.1 /zm and b' = 22 G/cm. 
The magnetic and optical trap were switched off simulta- 
neously and the atomic density distribution was recorded 
by absorption imaging along z after 18 ms of time-of- 
fiight (TOF) expansion. The atom number increases in 
steps of A'2 = 2 X 10'* every 8 kHz, in agreement with the 
expected Ai^2- We see three clear plateaus correspond- 
ing to 0, 1 and 2 sites spared from evaporation. For 
the shorter lattice period the frequency splitting is com- 
parable to the chemical potential. This results in some 
rounding off of the steps, but the plateaus remain visible. 

In the first set of experiments, we have characterized 
the free expansion of a single quasi-2D BEC [s^l- The 
clouds were released from the 5.1 /im period lattice and 
imaged after up to 18 ms of TOF. We extract the axial (/) 
and the radial (w) rms size of the cloud from gaussian fits 
to the density distribution. As might be expected, the 
observed expansion is predominantly one-dimensional. 
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Figure 2: Phase defects in quasi-2D condensates. Interference 
of four (a-c) and seven (d) independent BECs is observed 
12 ms after release from the 2.7 /im period lattice. Dislocation 
lines in the interference patterns (b-d) reveal the presence of 
phase defects in quasi-2D condensates. 

along the axial direction. For short expansion times, t < 
3 ms, the apparent axial size is limited by our imaging res- 
olution, but for longer times it follows the linear scaling 
I = vt, with V = 2.7mm/s. This value is comparable to 
the calculated velocity in the harmonic oscillator ground 
state along z, Vg — huj2 / {2m) — 2.1 mm/s, where m is 
the atomic mass. We find that the radial expansion can 
be described by the empirical law w — woy^l + {t/td)^, 
with wq = 4.4 /im and Iq — 5.7 ms. The same law de- 
scribes the radial expansion of a cigar-shaped 3D conden- 
sate, with ^0 = 2^ 1 "where lo^ — 1.5 ms for our trap. 
The radial expansion of our 2D gas is slower by a factor 
of « 4 compared to the 3D case, because the fast axial 
expansion results in an almost sudden (w^* = 76 /xs) de- 
crease of the atomic density, and only a small fraction of 
the interaction energy is converted into radial velocity. 

In the second set of experiments, we have studied inter- 
ference of independent quasi-2D BECs. Between two and 
eight clouds were released from the 2.7 /xm period lattice 
and allowed to expand and overlap [s^l ■ The resulting in- 
terference patterns were recorded by absorption imaging 
along the radial direction y. Due to the finite imaging 
resolution we observe only the first harmonic of the in- 
terference pattern with period ht/(mdi). Each image is 
thus the incoherent sum of the pairwise interferences of 
nearest-neighbor condensates. 

Interference of equally spaced, independent BECs pro- 
duces straight interference fringes (Fig. I3a}) as long as 
each BEC has a spatially uniform phase [24 . Issf . The 
main result of this paper is the observation of topo- 
logically different patterns, which reveal the presence of 
phase defects in quasi-2D condensates. Striking examples 
are "zipper" patterns (Fig.Efb)), where the fringe phase 
changes abruptly by tt across a dislocation line parallel 
to z. On both sides of the dislocation, the fringe con- 
trast is as high as in Fig. \^&) . We also observe "comb" 
patterns (Fig.EJc)), which show a dislocation with high 
fringe contrast on one side of the line, and vanishing on 
the other. Finally, we sometimes see "braid" structures 
with two dislocation lines (Fig.j^Jd)). Single dislocations 
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(zippers and combs) are clearly visible in about 15% of 
200 experiments with four interfering clouds |34l| . To ver- 
ify that the occurrence of defects is an equilibrium prop- 
erty of the system, we have checked that dislocations are 
still observed when holding the clouds in the lattice for 
500 ms after ramping down the gradient b'. 

The simplest phase configuration which can produce 
a sharp dislocation line is a s ingl e vortex in one of the 
condensates (see also [sE HE HE 113 )■ ^^^^^ 
of two interfering BECs, one can show that a centered 
vortex always leads to a zipper pattern (see a simula- 
tion of the expected pattern in Fig. |31 (a)). The zipper 
is indeed the only type of dislocation we clearly observe 
with two clouds. When more than two BECs interfere, 
the presence of a single vortex can result both in a zip- 
per and in a comb pattern, depending on the phases of 
the other condensates. In Fig.|21(b) we show a numerical 
simulation with four BECs leading to a comb. Increasing 
the number of interfering BECs enhances the probability 
that some of them contain defects ^3|> but the inter- 
pretation of images also becomes increasingly difficult. 
Further, for a large number of clouds, a single defect will 
not produce a clear dislocation line in the first harmonic 
of the interference pattern, because it affects only the in- 
terference with the two neighboring BECs. Already with 
four clouds, only half of 100 simulations with a vortex 
show clear zipper- or comb- type dislocations. The other 
half shows weaker dislocations which are not easily dis- 
tinguishable from straight interference fringes. 

Despite the agreement between simulations involving 
a vortex and the observed patterns, we point out that 
it is in general not possible to unambiguously deduce 
the underlying phase configuration from an interference 
image. For example, a dislocation line could also come 
from a dark soliton, where the phase of one of the BECs 
changes by tt across a line parallel to the imaging axis. In 
future experiments simultaneous imaging along a second 
radial direction could allow us to discriminate between 
different possible phase structures leading to the observed 
interference patterns. 

So far we could not observe a clear signature of vortices 
in images of single condensates taken along the axial di- 
rection z. We suspect that this is difficult because of the 
expansion properties of a 2D gas. Rotating 3D BECs, 
in which vortices are readily detected after TOF j ex- 
pand mostly radially, while our clouds expand mostly axi- 
ally. Therefore, any small misalignment with the imaging 
axis will significantly reduce the contrast. Interferometric 
detection along a radial direction offers a fundamentally 
superior signal, because a localized defect affects the ap- 
pearance of the whole image. 

It is important to assess the temperature of the clouds 
in which the observed phase defects appear. Precise ther- 
mometry at the end of the experimental cycle is diffi- 
cult, because the thermal cloud is very dilute. However, 
we can estimate lower and upper bounds for the tem- 
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Figure 3: Examples of numerical simulations of two (a) and 
four (b) interfering condensates. In both cases one randomly 
chosen BEC has a phase factor e**^ corresponding to a centered 
vortex, and the others have randomly chosen uniform phases. 
For simplicity, we model the clouds as gaussian wave packets 
and neglect interactions during the expansion. The images 
are convolved with a gaussian of 4 ^m rms width to simulate 
the finite imaging resolution. 

perature. During the 500 ms selective evaporation rou- 
tine, the clouds are heated due to three-body recombi- 
nation, and the only constant source of cooling is the 
finite lattice depth; atoms with an energy larger than 
Vq are accelerated away by the magnetic field gradient 
(Fig. ^b)). Assuming the largest realistic evaporation 
parameter rj = Vo/{kT) — 10, we get a lower bound for 
the temperature Tmin ~ 250 uK. To get an upper bound 
we note that at the beginning of the experiment the con- 
densed fraction is certainly above 50%. During the ex- 
perimental cycle the number of condensed atoms in the 
remaining sites drops by a factor of « 2. This means 
that, even if we neglect losses in the total atom number, 
the final condensed fraction cannot be less than 25%. 
Using the measured number of condensed atoms and in- 
tegrating the Bose distribution over the density of states 
in the lattice, we get Tmax ~ 500 nK. Since the num- 
ber of thermal atoms is different at Tmin and Tmax, the 
two bounds correspond to different condensation temper- 
atures Tc, and the estimated temperature range is more 
clearly expressed as 0.7 < T/Tc < 0.9. In this temper- 
ature range kT > /ia;i,2, so the thermal clouds are not 
fully in the 2D regime. 

The fact that the clouds are close to Tc is probably 
essential for the understanding of our observations, and 
a systematic temperature study will be the subject of 
future work. The probability for a thermal excitation 
of the system into a vortex state is cx e^^^*^^, where 
F ^ E — TS is the free energy associated with the exci- 
tation, E the energy, and S the entropy. Here we es- 
timate the conditions for F/kT to be of order unity. 
For a vortex in the center of the condensate, E ~ 
N [{hw^f/ij] ln(i?/0 [13, where N is the BEC atom 
number, R the size of the condensate, and ^ = h/ ^2ra^, 
the size of the vortex core. Equivalently, E/{kT) ^ 
(1/2) npA^ ln(i?/^), where no is the peak 2D atom den- 
sity and A is the thermal wavelength h j \/2wmkT . The 
number of distinguishable positions for a straight vortex 
of size ^ in a region of size i? is ^ i?^/^^, and the as- 
sociated entropy is Sjk ^ 21n(i?/^). In this estimate 
F/kT cx (noA^ — 4) vanishes for ripA^ — 4. In our ex- 
periment noA^ ~ 10 — 20 is a few times higher than this 
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value. However already this crude agreement suggests 
that the thermal excitation of vortices might be possible 
in our system. 

Thermal excitation of a tightly bound vortex- 
antivortex pair |l7l | is more likely than that of a sin- 
gle vortex. In that case the entropy is comparable and 
the energy is typically lower by the logarithmic factor 
ln(i?/^), in our case ^ 4. These pairs are difficult to 
detect with our interferometric scheme since they create 
only small phase slips in the fringe pattern. However they 
can play a significant role by screening the velocity field 
of a single vortex, thus lowering its energy and making 
its appearance more likely Q . 

The fact that ln(i?/^) is not large compared to 1 un- 
derlines the mesoscopic nature of our system. In a ho- 
mogeneous 2D system with i? — s- cx), both the energy and 
the entropy of a free vortex diverge as ln(i?), and the two 
contributions to the free energy cancel at the KT transi- 
tion temperature Tkt- Below Tkt only vortex- antivortex 
pairs are present, while above Tkt a large density of 
free vortices appears and suppresses superfluidity. In our 
case, we expect this phase transition to be replaced by 
a gradual increase of the average number of free vortices 
with temperature. For F ~ kT, the vortex number can 
show large fluctuations and two condensates produced 
under identical experimental conditions can have quali- 
tatively different wave functions. 

In conclusion, by selectively addressing individual sites 
of a ID lattice, we have produced both a single and sev- 
eral equally spaced quasi-2D BECs. We have character- 
ized the free expansion of a single BEC, and have inter- 
fcrometrically observed clear evidence for the presence 
of phase defects in about 10% of eight hundred conden- 
sates. While our observations can be explained by the 
presence of thermally excited vortices in the system, this 
does not exclude other scenarios and we hope that our 
experiments will stimulate further theoretical work. 
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